Phase ordering on small-world networks with nearest-neighbor edges 
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We investigate global phase coherence in a system of coupled oscillators on a small-world net- 
works constructed from a ring with nearest-neighbor edges. The effects of both thermal noise and 
quenched randomness on phase ordering are examined and compared with the global coherence in 
the corresponding XY model without quenched randomness. It is found that in the appropriate 
regime phase ordering emerges at finite temperatures, even for a tiny fraction of shortcuts. Nature 
of the phase transition is also discussed. 
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It is now known that diverse systems in nature may 
have the same topological structure as the small-world 
networks, first modeled by Watts and Strogatz (WS) [Q. 
A small-world network is usually constructed from a lo- 
cally connected regular network with given interaction 
range k, where some of the edges are randomly "rewired" , 
creating long-range "shortcuts". The small- world ef- 
fect in such a model refers to a high degree of clus- 
tering as well as the small characteristic path length 
i ~ logiV/log(2fc— 1), which is defined to be the aver- 
age of the shortest distance between two vertices in the 
network of size N (i.e., N vertices) [|]-|3). So far, a major 
of studies have mostly paid attention to the geometrical 
properties described by the above two quantities. On the 
other hand, some recent studies have considered dynam- 
ical systems put on small- world networks [^|-|6j, to find, 
e.g., faster propagation of information and better com- 
putational power. Such studies of dynamical systems, 
which apparently have wide applicability in physics, com- 
putational science, and biological science, have usually 
been performed on small-world networks with the inter- 
action range k > 2. The original WS model is poorly 
defined for k = 1 since the finite probability of isolated 
vertices induced by rewiring of the connections between 
the vertices yields divergence of the characteristic path 
length of the system. This difficulty may be circumvented 
by modifying the network construction method in such 
a way that shortcuts are added without removing local 
edges . This construction method leads to an increase 
of the number of total edges in the system; this is in 
contrast with the WS model where the total number of 
edges is conserved. 

In this paper we study a set of nonlinear coupled os- 
cillators on small-world networks constructed from rings 
with nearest neighbor edges (fc = 1), with attention to 
the emergence of long-range phase ordering as the num- 
ber of shortcuts is increased. In particular, the effects 
of thermal noise together with quenched randomness on 
phase ordering are explored and the aspects of the global 
coherence are compared with those in the corresponding 
XY model without quenched randomness. 



The small-world network considered in this paper 
is constructed in the following way: First, a one- 
dimensional regular network with nearest neighbor con- 
nections between N vertices is constructed, with the pe- 
riodic boundary conditions imposed. Then shortcuts are 
added with the probability P, between randomly chosen 
pairs of vertices. While local connections remain intact, 
more than one edge between any two vertices as well as 
any edge connecting a vertex to itself are prohibited, al- 
though such multi-connections and self-connections have 
been shown not to change the qualitative behavior of the 
WS model 0. Throughout this paper the probability P 
is defined to be the ratio of the number of total long- 
range connections (shortcuts) to the number of the total 
local edges in the network. 

At each vertex of this small-world network is located 
an oscillator; an edge connecting two vertices represents 
coupling between the two oscillators at those two ver- 
tices. Describing the state of the «th oscillator located 
on vertex i by its phase (jji , we write the set of equations 
of motion governing the dynamics of the N oscillator sys- 
tem (i = 1,2, • • ■,N): 

ij)i(t) =uji- J ^ sin(0j - + 77i(i), (1) 

where the intrinsic frequency uji of the ith oscillator are 
quenched random variables distributed according to, e.g., 
the Gaussian distribution with variance a^. The set A^ 
denotes the neighborhood of vertex i, consisting of the 
vertices connected to vertex i (via either local edges or 
shortcuts) with the coupling strength J. The last term 
rji on the right-hand side represents thermal noise, i.e., 
independent white noise with zero mean and correlations 

where the noise level T (> 0) may be regarded as the 
temperature of the system with the Boltzmann constant 
ks set equal to unity. When all the oscillators are iden- 
tical with no quenched randomness ((T^ = 0), Eq. (|l|) 
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describes the system of classical XY spins, for which the 
Hamiltonian is given by 



H 



(2) 



For example, niay represent the angular direction of 
the two-dimensional spin or the phase of the supercon- 
ducting order parameter at vertex i. The investigation 
of the phase transition in the XY model on a WS small- 
world network has revealed the presence of long-range 
order at finite temperatures, even for a tiny fraction of 
shortcuts On the other hand, in the absence of the 
thermal noise (T = 0), Eq. (|^) describes the dynamics of 
a set of coupled oscillators, which has also been studied 
on a WS small-world network and found to exhibit col- 
lective synchronization effectively ^ . These two systems 
can be described by Eq. in the appropriate limits; in 
this respect Eq. (|l|) provides a very interesting model, 
which allows to explore both the dynamic synchroniza- 
tion transition driven by quenched randomness and the 
thermodynamic phase transition by thermal noise to- 
gether. When the set includes the whole N oscillators, 
Eq. describes the dynamics of the fully-connected os- 
cillators with the normalized coupling strength J/N. The 
phase synchronization in such a globally coupled system 
has been extensively studied and the increase of the 
critical coupling strength due to the thermal noise has 
been found 

We here investigate the effects of quenched randomness 
together with thermal noise on the small-world network 
constructed as above. Collective behavior of the coupled 
oscillators on such a small-world network is conveniently 
described by the order parameter 
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where (• • •) represents the thermal average with respect 
to the thermal noise. For a small-world network, it is also 
necessary to take the average over different network real- 
izations, denoted by [•••]. To compute the order param- 
eter, we have integrated numerically the set of equations 
of motion given by Eq. (|^) through the use of Heun's 
method with the discrete time step At — 0.05. 
Typically, while the equations have been integrated for 
Nt = 4x 10'^ time steps, the data from the first Nt/2 steps 
have been discarded in measuring quantities of interest. 
Both At and Nt have been varied to verify that the mea- 
sured quantities are precise enough and the networks of 
various sizes, up to iV = 1600, have been considered. For 
each network size, we have performed one hundred inde- 
pendent runs with different configurations of the intrinsic 
frequencies as well as different network realizations, over 
which averages have been taken. 

Figure |l| displays the obtained order parameter m ver- 
sus the rescaled temperature T/ J at various values of the 



fraction P on a network of size N = 800 (a) in the ab- 
sence of the quenched randomness (ct^ = 0) and (b) in the 
presence of a finite amount of randomness (cr^ = 0.05). 
Thus Fig. |^(a) exhibits the behavior of the magnetiza- 
tion with temperature in the corresponding XY model. 
In the high-temperature limit (T/J — > oo), phases of 
the oscillators are distributed uniformly on the interval 
[0, 27r), leading to the absence of macroscopic coherence 
characterized by m = 0(A^^^/^). In the opposite low- 
temperature limit {T/J — > 0), all the oscillators become 
synchronized for P 7^ to give m = 1, regardless of the 
detailed structure of the network. On the other hand. 
Fig. ^b) indicates that the presence of randomness tends 
to suppress synchronization and macroscopic coherence 
is not achieved even at low temperatures for small but 
finite values of P. Indeed at low temperatures, the value 
of the order parameter for P < 0.2 in Fig. [|(b) as well as 
that for P = in Fig. |^(a) has been observed to diminish 
as the system size N gets larger, suggesting the absence 
of coherence (m = 0) in the thermodynamic limit. For 
P ^ 0.3 in Fig. |l](b), in contrast, the order parameter at 
low temperatures tends to increase with the system size, 
indicating the emergence of synchronization. 

It is thus manifested in Fig. |l| that phase ordering ex- 
hibits strong dependence on the fraction P. Regardless 
of the quenched randomness, in particular, the phases do 
not order in the absence of shortcuts (P = 0), which is 
consistent with the known result in one dimension [^ . 
When some fraction of the shortcuts comes into the sys- 
tem, on the other hand, the dynamics of the system 
changes dramatically, giving rise phase ordering. Here 
the critical value of the fraction beyond which the 
ordering emerges grows as the randomness is increased. 
Namely, larger fractions of the long-ranged connections 
are needed for ordering in the system with stronger 
quenched randomness. For instance. Fig. |l](a) shows the 
emergence of ordering even for P = 0.1 at low temper- 
atures, which implies Pc ^ 0.1. On the other hand, or- 
dering in Fig. 0(b) is observed for P > 0.3 (but not for 
P ^ 0.2 even at low temperatures), apparently indicat- 
ing that 0.2 Pc ^ 0.3. In this case of finite randomness 
(ct^ = 0.05), extensive simulations yield Pc ~ 0.23 in the 
low-temperature limit. Note that the critical fraction Pc 
in general increases with the variance cr^ and the tem- 
perature T, reflecting that quenched randomness as well 
as thermal noise tends to suppress synchronization. 

To explore the nature of the phase transition and to 
determine accurately the critical temperature Tc, we ex- 
amine the behavior of the order parameter by means of 
the finite-size scaling analysis. Recalling the critical be- 
havior of the order parameter in the thermodynamic limit 



m ^ (Tc - Tf 
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with the critical exponent (3, we expect the finite-size 
scaling form O 



(5) 
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with an appropriate scaling function F, where the criti- 
cal exponent D describes the divergence of the correlation 
volume at 
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\T-Tr. 



(6) 



According to this scaling, the plot of mN^^'^ versus T 
should give a unique crossing point at Tc. After /3/i/ and 
Tc are determined from the plot of mN^^'^ versus T, one 
then use 
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to obtain the value of (1 — (3)/v\ this, combined with the 
known value of finally gives the values of /3 and D. 

Figure ^ shows the determination of Tc on the small- 
world network with k — 1 and P — 0.5 in the absence 
of the quenched randomness (cP' = 0). Varying the 
value of we find that P/v k, 0.25 gives the well- 

defined crossing point at Tc ~ 1.02 (in units of J). In 
the inset of Fig. g, the least-square fit to Eq. (0) gives 
(1 - w 0.24, which, combined with j3/v « 0.25, 

yields (3 « 0.51 and v k, 2.04. These results demonstrate 
that the small-world network constructed from rings with 
nearest neighbor edges is also a mean-field system, which 
has the values /3 = 1/2 and v — 1/2 Further, the ob- 
tained value of D close to 2 indicates that the upper criti- 
cal dimension of the phase ordering transition is four [p^ , 
the same as that on the usual small- world network of the 
WStype§. 

Similarly, the scaling form in the presence of finite ran- 
domness is plotted in Fig.l3| displaying the determination 
of Tc and exponents for a — 0.01. Via the same analysis 
as in Fig. ^, we obtain Tc « 0.96 (again in units of J), 
which is smaller than the value in the system without 
quenched randomness. Thus the quenched randomness 
tends to suppress phase ordering, lowering the critical 
temperature. The corresponding values of the exponents 
/3 w 0.49 and v k, 1.96 are essentially the same as those in 
the absence of randomness (cr^ = 0). This concludes that 
the coupled oscillator system on a small-world network 
with the number of connections given by 0{N) displays a 
mean-field synchronization transition even for the short- 
est local- interaction range fc = 1, similarly to the system 
on a globally connected network with the much larger 
number of connections 0{N^). 

Figure || displays the phase diagram in the space of the 
rescaled standard deviation cr/J, the rescaled tempera- 
ture T/ J, and the logarithmic fraction — log P. The data 
for the boundary surface separating ordered (O) and dis- 
ordered (D) phase have been obtained from the finite-size 
scaling analysis, with the temperature T varied for given 
variance and the fraction P. In the absence of both 
quenched randomness and thermal noise (cr^ = T = 0), 
the frequencies of the oscillators are already ordered even 
without any long-range connections, and the phases can 
be ordered through the coupling between the oscillators. 
Therefore we have the critical fraction Pc = in this 



case, which is implied in the first and the second insets 
of Fig. IJ. The third inset displays the boundary on the 
plane of P = 1, which exhibits the interplay of the two 
types of the randomness (thermal noise and quenched 
randomness) in suppressing the phase ordering. 

In summary, we have examined phase ordering in the 
system of coupled oscillators on small-world networks 
constructed from rings with nearest-neighbor edges. The 
phase ordering has been found to depend strongly upon 
the fraction P of the long-range connections to the short- 
range local ones. In particular, such ordering, which 
is absent at P = for nonzero quenched randomness 
(ct^ ^ 0) , emerges at finite temperatures as P is increased 
from zero. The critical value Pc, below which ordering 
does not set in, apparently vanishes in the limits T — > 
and 0, increasing from zero with the randomness. 

The phase boundary between the ordered state and the 
disordered one, obtained from the finite-size scaling anal- 
ysis, has been found to be of the mean-field type. Finally, 
we note that the disordered state in Fig. || represents 
phase disorder, thus the possibility of a different type of 
order such as frequency order may not be excluded, the 
investigation of which is left for further study. 
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FIG. 2. Finite-size scaling of the phase order parame- 
ter [see Eq. ^] on the small-world network with k = 1 and 
P = 0.5, in the absence of quenched randomness (cr^ — 0). 
There exists a unique crossing point at Tc ~ 1.02. Inset gives 
the slope (1 — j3)/u « 0.24, which, combined with P/v ^ 0.25 
found in the main panel, results in /3 « 0.51 and v « 2.04. 
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FIG. 1. Phase order parameter m versus the rescaled tem- 
perature T/J on the small-world network with the local in- 
teraction range fc = 1 and size A'^ = 800 for the variance (a) 
= 0.0; (b) cr^ = 0.05. The error bars estimated by the stan- 
dard deviation have approximately the sizes of the symbols 
and the lines are merely guides to eyes. 
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FIG. 3. Finite-size scaling of the phase order parame- 
ter on the small- world network with fc = 1 and P = 0.5 in 
the presence of a finite amount of randomness (cr = 0.01). 
There is given a unique crossing point at Tc ~ 0.96, which 
is apparently smaller than T^ ~ 1.02 for — 0. Inset: The 
slope (1 — /3)/i/ « 0.26 is obtained, which, combined with 
P/f ^ 0.25, yields (3 « 0.49 and ^ 1.96. 
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FIG. 4. Phase diagram of the osciUator system on a 
smaU- world network. The data points on the phase boundary 
have been obtained from the finite-size scaling of the order pa- 
rameter in Eq. (^. O and D represent the ordered state and 
the disordered one, respectively. Insets: Phase boundaries on 
the planes T — 0, a — 0, and P = 1, respectively. 
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